Soluble Models of Strongly Interacting Ultracold Gas Mixtures in Tight Waveguides 
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A Fermi-Bose mapping method is used to determine the exact ground states of several models 
of mixtures of strongly interacting ultracold gases in tight waveguides, which are generalizations of 
the Tonks- Girardeau (TG) gas (ID Bose gas with point hard cores) and fermionic Tonks-Girardeau 
(FTG) gas (ID spin-aligned Fermi gas with infinitely strong zero-range attractions). We detail 
the case of a Bose-Fermi mixture with TG boson-boson (BB) and boson-fermion (BF) interactions. 
Exact results are given for density profiles in a harmonic trap, single-particle density matrices, mo- 
mentum distibutions, and density-density correlations. Since the ground state is highly degenerate, 
we analyze the splitting of the ground manifold for large but finite BB and BF repulsions. 

PACS numbers: 03.75.-b,03.75.Mn 



Due to the rapidly increasing sophistication of exper- 
imental techniques for probing ultracold gases, theoret- 
ical emphasis has shifted from effective field approaches 
to more refined methods capable of dealing with correla- 
tions, and experiments measuring such correlations have 
been carried out. In ultracold gases confined in de Broglie 
waveguides with transverse trapping so tight that the 
atomic dynamics is essentially one-dimensional with 
confinement-induced resonances l2| allowing Feshbach 
resonance tuning [3( of the effective ID interactions to 
very large values, such correlations are greatly enhanced. 
This has allowed experimental verification \4\, |5[ of the 
fermionization of bosonic ultracold vapors in such geome- 
tries predicted by the Fermi-Bose (FB) mapping method, 
which was introduced in 1960 [6j and used to obtain the 
exact iV-particle ground and excited states of a ID gas 
of impenetrable point bosons, the Tonks-Girardeau (TG) 
gas. Apart from recent work of Imambekov and Demler 
7], theoretical work on ID mixtures has used approxima- 
tions which are not valid in the highly-correlated regime. 
Exact results in [7| , based on the Bethe ansatz method of 
Lieb and Liniger (LL) Q , are limited to the spatially uni- 
form case (no longitudinal trap potential). We present 
here exact solutions for several models of mixtures of 
strongly correlated ID systems including that of [7[, for 
both untrapped and harmonically trapped mixtures. 

Bose-Fermi mixture: We consider a ID mixture of 
Nb bosons with point hard-core boson-boson interac- 
tions and of Np noninteracting fermions with point hard- 
core boson-fermion interactions, and we assume equal 
masses = mp = m. This model could be realized 
by choosing a bosonic and fermionic isotope of a given 
alkali element (e.g. 6 Li - 7 Li). We indicate the boson and 
fermion coordinates by Xb — (a^is,--- ,xn b b), Xp — 
(xip, • • ■ , xn f f)- Our treatment includes the case where 
the mixture is subjected to an external potential (e.g. a 
harmonic trap v(x) — ^muj 2 x 2 ), provided that the po- 
tentials acting on the bosons and on the fermions are the 
same, i.e., vb(xjb) — v{xjb) a,ndvF(xep) — v{xip). The 



Schrodinger Hamiltonian is H = Hb ■ 
with 



-Hf + Hbb + Hbf 



Nb 



i=i 

Hp = J2[-(h 2 /2m)d 2 /dxj F + v(x lF )} ■ (1) 
i=i 

The impenetrable point BB and BF interactions can be 
formally represented as sums of LL delta function inter- 
actions |8J gBB^(xjB — %ib) and gBF${xjB — xif) with 
coupling constants gsB —> +°o and qbf —>■ +oo, but 
they are more conveniently represented by constraints 
that the many-body wavcfunction ^>(Xb,Xf) vanishes 
at all BB and BF collision points, after which Hbb and 
Hbf can be omitted from the Hamiltonian. 

The exact solution can be obtained as follows. Con- 
struct a "model wavefunction" $>m(Xb,Xp) which is a 
Slater determinant ofN= Ns+Np orthonormal orbitals 
u\(x), • • • , un(x) occupied by all N bosons and fermions, 
with all possible permutations of the atoms among the 
orbitals: ^ m = Yip £ {P) u i{Px\) • ■ • un(Pxn)- Here 
xi, ■ ■ ■ ,xn are x\b, ■ • • ,xn b b,xif, • • • , xn f f, the sum 
runs over all AH possible permutations of these variables 
including permutations exchanging bosons with fermions, 
and e(P) is the usual ±1 sign of the permutation. This 
wave function vanishes not only at the points Xjp = X£p 
required by fermionic antisymmetry, but also at the 
points XjB = X£b and XjB = x^p required by the BB 
and BF hard core constraints. Its improper symmetry 
under BB and BF exchange is repaired by a modified FB 
mapping. Define a mapping function A{Xb, Xp) by 
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i<j<e<N E 



N B N f 



Sgn(x jB—XiB 



i Y[ Y[sgn(x jB -x lF ) (2) 



where the sign function sgn(x) is +1 ( — 1) if x > 
(x < 0). Then the physical wavefunction is 
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^(X B ,X F ) = A(Xb,X f )^ m (X b ,X f ). * m is an ex- 
act many-body energy eigenstate if the orbitals u v are 
eigenfuctions of the single-particle Schrodinger Hamilto- 
nian — (h 2 /2m)d 2 /dx 2 + v(x) with eigenvalues e v . Since 
A is constant except for jumps at nodes of fifi it 
follows that the physical state \& is an energy eigen- 
state with eigenvalue J2 U €v ano - ^ s symmetric under 
permutations of bosons and antisymmetric under per- 
mutation of fermions. The ground state is a filled 
Fermi sea of the lowest N orbitals and excited states 
are generated by choosing higher orbitals. An impor- 
tant case for experiments is that of harmonic trap- 
ping, v(x) = ^muj 2 x 2 . Then the ground state is a 
straightforward generalization of that of the trapped TG 
gas Q. The orbitals are Hermite-Gaussians (j) n {x) — 
7r- 1 / 4 (2 n n\x osc )- 1 / 2 e-Q 2 / 2 H n (Q) with n = 0,1,- •• ,N, 
where Q — x/x osc and a; osc = y/h/mu}. ^om can be 
reduced to Bijl-Jastrow form by determinantal algebra 
as in the original solution 0] for the untrapped TG gas 
and the recent one [§] for the trapped TG gas, yielding 

*OM °c [I\i<j<e<N( x 3 ~ x *)] HjLi <Po( x j) where xi = 

%1B,--- ,%N B — X N B B,XN B + 1 = X\F,--- > %N = %N F F- 

By applying the mapping function A to \&oAf one finds 
for the ground state 



N b N F 

*0 = [ II \XjB- X tB I] [JI \XjB - XlF I] 
l<j<£<N B j=lt=l 

X [ II {x ]F - X lF )]<j) {xiB) ■ ■ ■ <Po{x NbB ) 

i<j<e<N F 

xcj) Q (x 1F ) ■ ■ ■ (f>o(x NF p) . (3) 

apart from normalization. The solution for the homoge- 
neous system of size L with periodic boundary conditions 
differs only by omission of the Gaussian factors </>o and 
by replacement of each coordinate difference x — x by 
sin[7r(a; — x')/L] [6]. The bosonic and fermionic reduced 
one-body density matrices are defined as P\b{x,x') = 
m B J^ (x,X' B ,X F )%(x',X B ,X F )dX' B dX F and 
p 1F (x, x') = yi F J ^ (XB,X F ,x)^* (XB,X' F ,x')dX B dX' F , 
where X' B = (x 2 b, • • ■ , x NbB ) and X' F = 
{x\ F ,--- ,xn f -i,f), and 91b, Vtp are normalization 
constants fixed by the conditions J p\B{x,x)dx = Nb 
and J p\ F {x,x)dx = N F . 

Density profiles: The bosonic and fermionic single- 
particle densities Pb{x) = p\b{x,x) and Pf{x) = 
p\ F {x, x) are the same as those of $om, since they de- 
pend only on |^o| 2 i and A 2 = 1. Since ^qm is com- 
pletely antisymmetric under permutations of all N par- 
ticles, ps and p F are both proportional to the density 
Ptg °f a harmonically trapped TG gas of N bosons: 
p B {x) = (N B /N)p TG (x) and p F (x) = {N F /N)p TG (x) 
where Q Ptg(x) = J2n=o \4>n{x)\ 2 . These exact results 
contrast strongly with local density approximation re- 
sults in the TG limit 0], which show BF phase separa- 
tion with ps more concentrated in the center, although 



their total density pb(x) + p F (x) agrees closely with ours. 
Note that the exact ground state is highly degenerate in 
the TG limit [l(| due to the fact that there is no re- 
quired symmetry under BF exchange. Our choice ([3|) is 
the most symmetrical, has hard-core cusps at all BB and 
BF collisions and antisymmetry nodes at all FF colli- 
sions, and it will be shown below that if the degeneracy 
is lifted by making the BB and BF interactions finite, 
then ([3]) corresponds to the lowest member of the split 
ground manifold, and hence our pb{x) and p F {x) are the 
TG limit of the true ground state partial densities for 
large but finite repulsions. Nevertheless, in the TG limit 
the density profiles are labile, since linear combinations 
of the degenerate ground basis [l0( have different density 
profiles. 

Density- density correlations and collective excitation 
spectrum: The BB, FF, and BF density-density corre- 
lation functions for this model depend only on \l/oM) and 
hence are all the same as that of an ideal Fermi gas or 
TG gas of N particles apart from normalization; see Eqs. 
(10) and (11) of [9]. This implies, e.g, that no composite 
fermions are found (as is the case for a BF mixture on a 
lattice for particular values of gsF and qbb H])- The 
spectrum of collective excitations is given by the poles of 
the dynamic structure factor, and hence coincides with 
that of an ideal Fermi gas. For the case of harmonic 
trapping the spectrum is given by integer multiples of 
the trapping frequency to for both in-phase and out-of- 
phase modes, in disagreement with 0]. 

Momentum distributions: Momentum distributions, 
Fourier transforms of the reduced one-body density 
matrices with respect to x — x', differ greatly be- 
tween bosons and fermions. Consider first the bosonic 
one-particle density matrix p\ B [x,x') and momentum 
distribution ns(fc). Using (O one finds that p\B is 
proportional to the one-body density matrix of a pure 
TG gas of N bosons, p\ B {x,x') = (Nb/N)p\tg{x,x'), 
and hence ns{k) = (NB/N)riTG(k)- Although 
Pytg is not known analytically, the asymptotics of 
nTc{k) are known for a homogeneous system of den- 
sity pre and are given by nTc(k) oc |fc| -1 / 2 for 
k <C pre UM an d by nrcik) oc |fc|~ 4 for k ^> pro 



[13| . Determination of p\ F (x,x') and n F (k) is much 
more difficult. We detail here the derivation for 
the homogeneous system. Comparing Eq. with 
of the pure TG gas Q we obtain pi F (x, x') — 
91f f dxi...dx N -i\^ Q MM-i\ 2 n i < ::j < NB sgn(xj 
x)sgn(xj - zOni^iV-i 4 ^^!^ ^ x)/L)sm(Tr(xj - 
x')/L) where ^o&r.iV-i is $om for N — 1 parti- 
cles. The determinantal expression $om,jv-i = 
det[e* fc ' ! %^ = i i ..jv-i ) where kj = -(N - 2)n/L, —(N - 
4)n/L, (N - 4)tt/L,(N - 2)n/L and some alge- 
bra reduce p\ F to the sum of a product of one- 
dimensional integrals, of which the first Nb involve 
b(t,x,x') = |sin(7r(i - x)/L)\\sm(-!r(t - x')/L)\ 
as in a TG gas, and the remaining N F contain 
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f(t,x,x') = sin(7r(£ — x)/L)sm(ir(t — x')/L) as in an 
ideal Fermi gas. Combining the dissimilar factors by the 
"phase trick" of @, 01 , one finds the final result 



d<t> 



2tt 



Pif(x,x') = a 
where a = (N F /NL)C{N B ,Nr) 

4 rL 



[u jtk (x,x',<j))] 



iNb<p det 

j,k=l,...N-l 

(4) 

and Uj : k(x, x' , cf>) = 
| fa A[e^sgn(t - x)sgn(t - x') + 1] sin(7r(t - 
x)/L) sin(7r(f - x')/i)e l27rt( ^ fc)/i . This can be ex- 
tended to a harmonic trap. 

Ground and excited states for large but finite BB and 
BF repulsion: Suppose now that the interactions are fi- 
nite LL interactions [8( gBB&ixjB —xib) and gBF$(xjB — 
xif). Define dimensionless coupling constants Xbb = 
"WV^t ^d Xbf = y/$=. The TG limit, to 
which our ground state corresponds, is Xbb +°° 
and Xbf ^ +oo. As previously noted this state is only 
one member of a degenerate ground manifold (loj . and 
we now examine how this degeneracy is split when Abb 
and Xbf are large but finite. The exact solutions of 
the two-body problem are known [lq and the energy 
eigenfunctions ^ v {xbf) of the relative motion in a har- 
monic well are expressible in terms of parabolic cylinder 
functions D v {^) [16[ where £ = xsF\p^f- The corre- 
sponding relative energy eigenvalues are E v = (v + ^)huj 
where v is a solution of the transcendental equation 
r(i(l — v))/T(— hv) = —Xbf, and reduce to the har- 
monic oscillator (HO) quantum numbers n — 0, 1, 2, • • ■ 
in the nonintcracting limit. Here we want the opposite 
limit Abf > 1. In the TG limit all the eigenfunctions 
are ^ti(xbf) = A(xbf)^um(xbf) where A is a map- 
ping function [10( and the model functions ^ n M are the 
odd HO functions u n with n = 1, 3, 5, • • • so as to vanish 
at contact. There are two possible mapping functions: 
A\ = sgn(x bf) and A2 = + 1, and the ground state 
is twofold degenerate, with even and odd eigenfunctions 
"L^ = sgn(xB_F)ui(xBF) and 9^ = A2U1 = u\(xbf) and 



energy E\ 



\tno. For large but finite Xbf the lowest 



even eigenfunction is — D Vl {\(\) where v\ is the so- 
lution of the above transcendental equation closest to 1. 
fyf does not vanish at xbf = and has an LL cusp 
[8| there due physically to the delta function interaction 
and mathematically to the absolute value in its argument 
171 ] . Using gamma function identities and Taylor expan- 
sions one finds v\ = 1 — + 0(Xg 2 F ). On the other 
hand, the lowest odd eigenfunction is the lowest odd 
unperturbed HO eigenfunction u\, since u\8(xbf) = 0. 
Hence the symmetric solution is the ground state 
and = i*i is the first excited state, with excitation 



energy AE = 



O(Xgp). More generally, all of 



the eigenfunctions in the order of increasing energy are 
u 1j -C^G^Qj m 2, • • • > alternately even and odd 
with 0, 1, 2, ■ ■ • nodes. This is in agreement with Sturm- 
Liouville theory and the theorem that the ground state 



of a Boltzmann or Bose system is real, nonnegative, and 
hence nodeless. 

For arbitrary Nb and Np and finite coupling constants 
the exact solutions are known in the untrapped case only 
for the special case Abb = Abf- However, the above 
theorem on ground state symmetry has a generalization 
according to which the ground state is nodeless except for 
Fermi antisymmetry nodes. It follows that the ground 
state vanishes only at FF contact points Xjs — Xif, and 
at all BB and BF collision points it has interaction cusps 
of LL form Q. By comparison with Eq. ©, one con- 
cludes that is the TG limit of the finite-interaction 
ground state; all other choices [lj| of the TG mapping A 
correspond to excited states. The manner in which the 
degeneracy of the TG ground manifold splits when Abb 
and Abf are made large but finite is related to symme- 
tries of all the different mappings [l(| yielding the differ- 
ent ground states degenerate in the TG-limit. Consider 
first a mapping A^ differing from ([2]) by omission of all 
of the factors sga(XjB— xipy The corresponding mapped 

state is = A^^qm where foM is the previously- 
defined fermionic model ground state. A^ restores the 
correct symmetry under BB permutations and introduces 
the necessary interaction cusps at BB collision points [j| , 
while leaving antisymmetric not only under FF per- 
mutations, but also locally antisymmetric under BF ex- 
changes, in the sense that changes sign as XjB — xgp 
passes through zero []J]. It follows that * 0) is an eigen- 
state not only of the TG Hamiltonian Abb — > 00 and 
Abf - * 00, but also of a Hamiltonian with TG BB in- 
teractions (i.e., Abb — * 00) but finite BF interactions, 
because the nodes at BF collision points "kill" the delta 
function interactions, and its energy is the same as that 
in the TG limit. It is not the ground state because of the 
previously-stated theorem, and in fact it is the top of a 
ladder of closely-spaced levels extending upward from the 
ground state, with splittings of order between adja- 
cent levels, obtained by restoring more and more collision 
sgn factors in the mapping function and hence collision 
cusps in the wave function, thus generating additional 
downward shifts of order -SSsL, 

7b f 

Other soluble models : In addition to the previously 
discussed Bose- Fermi mixture (which we denote as Model 
I) we can solve several other models with a similar ap- 
proach. Model II is a mixture of two noninteracting 
Fermi gases with interspecies point hard-core repulsions, 
and model III is a mixture of two Bose gases A and B 
with inter- and intraspecies hard-core point interactions. 
Model II could be realized, e.g., by choosing two hyper- 
fine levels of 40 K or 6 Li and by tuning to very large and 
repulsive the interspecies interactions by a confinement- 
induced Feshbach resonance and analogously model III 
could be realized by choosing two hyperfine levels of a 
bosonic alkali atom e.g., Na or Rb. Models I-III all start 
from the same model state ^m, so we only state final 
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results, again for harmonic trapping. For both models 

II and III the density profiles are proportional to those 
of a TG gas of N = Na + Nb bosons, and in model 

III the momentum distributions for both components are 
also proportional to those of a iV-particle TG gas, while 
for model II the ^4-component momentum distribution 
riyt(fc) is equal to np(k) for model I with Na fermions 
and Nb bosons, and similarly for 7ig(fc) with A and B in- 
terchanged. These models have degenerate ground states 
by the same mechanism as model I. 

Finally, we can consider binary mixtures with large, at- 
tractive, odd-wave interactions between the two species 
due to a p-wave Feshbach resonance, in the so-called 
fermionic Tonks-Girardeau (FTG) limit 0, [H, [I| 
wherein the fermionic wave function is nonzero with sign 
change at contact. Model IV is a mixture of two Bose 
gases A and B with no AA or BB interactions but an AB 
interaction of FTG form, model V is a mixture of ideal 
Bose gas B and FTG gas F with FTG BF interactions, 
and model VI is a mixture of two FTG gases A and B 
with FTG AB interactions. These models all have the 
same model state ^m, a mixture of two dissimilar ideal 
Bose gases, while the mapping function A has sgn fac- 
tors for all collisions where an FTG interaction is desired 
, [lii , l2(ij ]. The momentum distributions for the spa- 
tially uniform case on a ring have Lorentzian form as for 
the FTG gas [2l|, with the range of the density matrix 
of each component determined by the sum of the densi- 
ties of all components (either the same or different) with 
which it interacts. 

We have identified further models with attractive in- 
terspecies FTG interactions solvable by mapping from 
ideal gas mixtures to mixtures of an ideal Fermi gas or 
TG Bose gas and an ideal Bose gas or FTG Fermi gas, or 
a mixture of two dissimilar ideal Fermi gases or TG Bose 
gases, and we believe that this covers all two-component 
mixture models solvable by mapping functions every- 
where ±1, with jumps at contact points which convert 
Fermi antisymmetry nodes to Bose hard core cusps and 
Bose-Bose or dissimilar species contact points to infinite 
FTG attraction discontinuities. 

Experimental realizability and applications: Since the 
TG regime in pure ultracold gases has been reached 0, Q 
and p-wave resonances required for the FTG regime have 
been observed [22J , several of these strongly interacting 
mixture models should be experimentally realizable. The 
low-lying level structure we have found for Bose-Fermi 
mixtures might be exploited for quantum computation. 
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